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Abstract —In the present paper, a simplified two-dimensional 
magnetic recording (TDMR) channel model is proposed in 
order to capture the qualitative features of writing and read- 
back processes of TDMR systems. The proposed channel model 
incorporates the effects of both linear interference from adjacent 
bit-cells and signal-dependent noise due to irregular grain bound¬ 
aries between adjacent bit-cells. The simplicity of the proposed 
model enables us to derive the closed form of the conditional 
PDF representing the probabilistic nature of the channel. The 
conditional PDF is Gaussian distributed and is parameterized by 
a signal-dependent covariance matrix. Based on this conditional 
PDF, a Monte Carlo method for approximating the symmetric 
mutual information of this channel is developed. The symmetric 
mutual information is closely related to the areal density limit 
for TDMR systems. The numerical results suggest that we may 
need low-rate coding, e.g., 2/3 or 1/2, when the jitter-like noise 
becomes dominant. 

I. Introduction 

Continuing demand for high-density magnetic recording 
requires novel technology for realizing higher areal den¬ 
sity. The emergence of shingled writing promotes research 
into advanced signal processing because narrowing inter¬ 
track pitch causes severe inter-track interference and jitter¬ 
like noise around track edges. Under such circumstances, two- 
dimensional signal processing based on multiple reading heads 
is becoming a hot research topic in this field. For example, such 
a system requires linear interference cancellation in both the 
cross-track and down-track directions. 

In the near future, the size of bit-cells will likely shrink to 
be comparable to grains of a magnetizing material. As such, 
highly advanced two-dimensional signal processing will be 
required to handle strong linear interference from adjacent bit- 
cells and signal-dependent noise due to irregular grain bound¬ 
ary between adjacent cells. The concept of two-dimensional 
magnetic recording (TDMR) advocated by Wood et al. 
has inspired research in this field. Related research into, for 
example, signal processing 0, detection algorithms [ll, and 
simplified channel models 0 has been conducted for TDMR 
systems. 

One goal of the present paper is to develop a simple 
mathematical model for TDMR systems that is useful for 
performance evaluation of detection algorithms and for the 
design of two-dimensional codes. The simplified TDMR chan¬ 
nel model presented herein incorporates the effects of both 
linear interference from adjacent bit-cells and signal-dependent 
noise. Although the channel model is fairly simple, the model 
shows the complicated nature of signal-dependent channels. 


Observing the behavior of the channel model will help to 
clarify the characteristics of TDMR channels. The second goal 
of the present paper is to present a method for evaluating 
the symmetric mutual information. The evaluation of mutual 
information is not a trivial problem for a channel with signal- 
dependent noise, but the problem is worth pursuing because it 
is closely related to the areal density limit of TDMR systems. 

II. Two BIT-CELL MODEL 

Before introducing the simplified TDMR channel model, 
we will discuss a two-bit-cell model in this section. Although 
this channel model is simple, it provides insight for a TDMR 
channel and this model has a number of features that a TDMR 
channel model must possess. For example, this channel model 
can handle signal-dependent noise generated by the irregular 
boundary between two adjacent bit-cells. 

A number of grains of a magnetizing material are spread 
over the surface of the recording medium. In the writing phase, 
a writing head magnetizes each grain upward or downward, 
depending on the information to be written. A magnetization 
process of each grain is determined according to the strength 
distribution of magnetic flux emitted from a writing head and 
the magnetic sensitivity of the grain. The first approximation 
introduced in the present paper is that the head footprint of 
magnetic flux has an exactly rectangular shape. Namely, the 
magnetic held induced by a writing head has an effect on 
the rectangular area. This approximation helps to simplify the 
channel modeling process described below. This rectangular 
area is referred to herein as a bit-cell. Some of the grains 
lie between two adjacent bit-cells and cause irregular grain 
boundary between adjacent bit-cells. Since a read-back system 
has no knowledge on the positions of grains, such an irregular 
grain boundary incurs jitter-like signal-dependent noise that 
makes it difficult to achieve reliable detection in TDMR 
systems. 

A read-back signal is obtained by taking the convolution of 
a sensitivity function (i.e., 2D-impulse response) of a reading 
head and a magnetic held pattern caused by magnetized grains 
on the surface of a recording medium. A sample sequence 
of the read-back signals is sent to a signal processing and 
detection unit to estimate the written data. In a TDMR system, 
two adjacent bit-cells in a grid of bit-cells are closely arranged 
and a reading head with a wide (compared to the size of a bit¬ 
cell) footprint is often considered in a TDMR system. Two- 
dimensional linear interference consisting of both inter-symbol 
interference (ISI) and inter-track interference (ITI) arises in 
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Fig. 1. Two bit-cell model 


such a system and should be appropriately processed to provide 
high reliability of estimated data. 


A. Details 


The two-bit-cell model, which is an abstraction of a part of 
a TDMR system, is shown in Fig. [T] This model contains two 
adjacent bit-cells, and the information written to these bit-cells 
is denoted by xi,X 2 S {-1-1,—1}. Between two bit-cells, we 
assume that an irregular grain boundary exists. We also assume 
that two overlapping footprints (denoted by Head 1 and Head 
2 in Fig. [TJ of a reading head are exploited to yield the two 
read-back signals. 


In the two-bit-cell model, the two read-back signals yi and 
j /2 are given by 


yi = xi + [-+nj ]xi + 


1 


2/2 = X2 + [ - — Uj ] X2 + { - + rij ] xi + n2- (1) 


1 


1 


— n, ] X2 + ni 


Note that the read-back signals are assumed to be appropriately 
scaled so as not to introduce additional scale parameters. The 
symbols ni and 712 are independent Gaussian random variables 
with mean 0 and variance that represent system noises. The 
symbol Uj is an independent Gaussian random variable with 
mean 0 and variance that corresponds to the irregular grain 
boundary between two bit-cells. In other words, we introduce 
a stochastic model for expressing the effects of an irregular 
grain boundary. 


The following is an implicit assumption to derive the above 
model of the read-back signals. Let f{a,b) be the vertical 
component of magnetic flux induced by the media at point 
{a,b) G K. For simplicity, it is assumed that f{a,b) = Xi G 
{-fl, —1} if (a, 5) is in the magnetized area corresponding to 
bit-cell 1; otherwise, f{a,b) = X 2 G {-fl,—1} is assumed. 
Under this setting, the read-back signal yi is obtained by 

yi= f ficL,b) da db. (2) 

J Head i 

A read-back signal can be considered as an output from a 
concatenated system of a 2-D low-pass filter with a rectangular 
window (i.e., head footprint) followed by a 2-D sampler. 
Although this assumption is based on an ideal case (such a 


2-D low-pass filter is not realizable), we expect the model to 
reflect the qualitative nature of an actual TDMR system. 

If there is no system noise or jitter-like noise due to the 
irregular grain boundary between two bit-cells, the following 
read-back signals can be obtained: 

3 1 3 1 

2/1 = -X1 + -X2, 2/2 =2^2+2^1- 

Some of the grains located in the middle of the bit-cell 
boundary have the same polarity as the value of xi, whereas 
others have the polarity same as the value of X 2 - The balance 
of contributions of such grains are modeled by a Gaussian 
random variable rij in the proposed model. We adopt a 
Gaussian model for representing this balance for two reasons: 
1) as described later herein, a Gaussian model leads to a 
mathematically tractable probability density function (PDF) of 
received signals, and 2) since the balance of areas Uj consists 
of contributions from several grains (i.e., addition of several 
random variables), it is natural to model the balance of areas 
rij by a Gaussian random variable. 

We can rewrite the model ([T]i in vector form, as follows: 



The first term on the right-hand side of the above equation 
represents a linear interference that occurs on a:i and X 2 , and 
the second term on the right-hand side is the system noise. 
The third term on the right-hand side indicates the jitter-like 
noise due to the irregular grain boundary. Note that the last 
term on the right-hand side disappears if xi = X 2 - Namely, 
when two bit-cells have the same polarity, an irregular grain 
boundary has no harmful effect regarding detection. The last 
term on the right-hand side, representing the signal-dependent 
noise, introduces a certain difficulty in designing detection 
algorithms, codes, and evaluating the capacity of the channel. 

Figure |2] indicates signal constellations regarding the two- 
bit-cell model. The transmitted symbols {xi,X 2 ) are depicted 
in Fig. 12] (left), and the received symbols ( 2 / 1 , 2 / 2 ) in a noiseless 
case are indicated as black dots in Fig. |2] (right). Due to the 
linear interference, we can see that the received symbols are 
located at the corner points of a diamond-shaped region. The 
dotted circles around the black dot represent the contours of 
the PDF PY{yi,y 2 ){Y = ( 2 / 1 , 2 / 2 ))- The exact form of the PDF 
will be discussed in the next subsection. Note that the shape 
of the contours around the dots depends on the signals. For 
example, the noiseless received signals ( 2 / 1 , 2 / 2 ) = (+2,+2) 
correspond to the transmitted signals {xi,X 2 ) = (+l,-fl). In 
this case, no jitter-like noise appears in the detection phase, 
and only system noise ( 77 , 1 , 77 , 2 ) disturbs the received signals. 
On the other hand, when {xi,X 2 ) = (-fl,—1), the irregular 
grain boundary is as depicted in Fig. |2] (right). 


B. Conditional probability density function of received signals 

One of the advantages of the two-bit-cell model is its math¬ 
ematical tractability. It can be shown that the conditional PDF 
of the received symbols are Gaussian distributed (details are 
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Fig. 2. Signal constellations (left: transmitted symbols {xi^X 2 ), right: 
received symbols (yi,y 2 ) in a noiseless case) 


discussed in the next section). The conditional PDF FY|x(y|x) 
(x = {xi,X 2 )'^,y= ( 2 / 1 , 2 / 2 )^) is given by 

Pv,x(y|x) = „ exp (-(y - Ax)^5(x)-i(y - Ax)) 


27r\/|S'(x)| 


where the matrix A is the interference coefficient matrix 
defined as 

I f 3 1 
1 3 


A = 


(4) 


The covariance matrix S (x) is given by 
S{x.) = f + (j 1 a']{xi-X2Y 


cf'^Axi - X2Y 


CTjixi - X2)^ + ol 


■ (5) 


for a channel model containing a number of bit-cells. Another 
observation we can make from this figure is that we may 
need a non-linear detector instead of a linear detector in 
order to approximate the exact ML detector with reasonable 
computational complexity. 



Fig. 3. ML decision region (horizontal and vertical axes represent yi(—6 < 
yi < -1-6) and 1/2 (—6 < y 2 < -1-6), respectively. The red dot indicates the 
received points in the noiseless case.) 


III. Simplified TDMR channel model 

In this section, we generalize the two-bit-cell model dis¬ 
cussed in the previous section into a more general two- 
dimensional channel model that is referred to as a simplified 
TDMR channel model. 


From FV|x(y|x), the noise vector y — Ax is Gaussian 
distributed, and the covariance matrix S'(x) is dependent on 
the transmitted signal x. This observation matches the contours 
shown in Fig. |2] (right). 

C. Maximum likelihood detection 

We have had the exact conditional PDF for the two-bit-cell 
model. The conditional PDF leads to the explicit formula for 
maximum likelihood detection for this channel. At a glance, 
the two-bit-cell model appears to be a simple model, but this 
model possesses properties that are not so intuitive. In this 
subsection, we will observe such properties in terms of signal 
detection. 

Let us define a distance measure Z3(y|x) as 

D{y\^) = (y - Ax)^5(x)-i(y - Ax) + log |5(x)|, (6) 

which is derived by taking the negative logarithm to the 
conditional PDF of the channel. Using this distance measure, 
the maximum likelihood (ML) detection rule can be written as 

X = arg^g{+i__i |2 minD(y|x). (7) 

Figure [3] shows the decision regions of the two-bit-cell 
model for this ML detection rule in the ( 221 , 2 / 2 ) plane. From 
left to right, the cases of standard deviations cjj = 0.01 
and 0.25 are depicted. Unlike a conventional linear inter¬ 
ference channel, we can see that two decision regions for 
( 221 , 2 / 2 ) = (+1,+1), (“1, “1) become closed sets as the 
standard deviation tjj increases. The peculiar shape of these 
decision regions suggests the difficulty of a direction problem 


In the following discussion, we will assume that bit-cells 
are arranged on grid points in a rectangular area, as shown in 
Fig. a The area of interest is called the writable region that 
consists of I rows and J columns. Each bit-cell has its index 
i{i £ [l,n],n = IJ), and the bit-cell with index i is denoted 
by Bi. As in the case of the two-bit-cell model, each bit-cell 
can take the value -|-1 or —1. The notation [a, b] represents the 
set of consecutive integers from a to b. 


J cells 




i 

1 

2 

3 

4 

5 


6 

7 

8 

9 

10 


11 

12 

13 

14 

15 


16 

17 

18 

19 

20 


21 

22 

23 

24 

25 


□ -1 
□ +i 


Fig. 4. Writable area and bit-cells {n = 25) 

The set of indices of bit-cells adjacent to the bit-cell Bi is 
denoted by N{i){i £ [1, n]). For example, in the case of Fig. a 
we have iv(13) = {8, 12, id, 18}. The boundary between two 
adjacent bit-cells is referred to as an edge. The edge between 
bit-cells Bi and Bj is represented by Ei^j, or equivalently Ej^i. 




















































The head footprint Hi corresponding to bit-cell Bi{i G [l,n]) 
provides a received symbol yi by reading a certain area around 
Bi. Here, we assume that the head footprint has a crisscross 
shape, as shown in Fig. |5] 
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Fig. 6. Expressing irregular grain boundary with Gaussian random variable 


Fig. 5. Head footprint Hi (gray area corresponds to the head footprint Hi) coefficients and q(x) denotes the noise term 


In the simplified TDMR channel model, the received signal 
yi{i G corresponding to the head footprint Hi is 

modeled by 

yi = axi + P ^ Xj + ^ \Xi- Xj\zij + Zi. ( 8 ) 

jGN{i) jGN{i) 

The real constants a and /3 are scale parameters for the 
intended signal and the linear interference coming from ad¬ 
jacent cells, respectively. The noise term Zi is a Gaussian i.i.d. 
random variable with mean 0 and variance and represents 
a system noise. An i.i.d. random variable Zij is identical 
to Zj^i and represents jitter-like noise due to the irregular 
grain boundary around edge Eij. The random variable Zij 
follows the one-dimensional Gaussian PDF with mean 0 and 
variance Of course, it is easy to provide more flexibility 
by introducing more parameters. For example, the coefficients 
of linear interference from the adjacent cells can be direction 
dependent. However, we will avoid introducing too many 
parameters so as not to complicate the following discussion. 
Once the following argument is grasped, modification of the 
model would be straightforward. 

As in the case of the two-bit-cell model, jitter-like noise 
occurs around an edge between two bit-cells having distinct 
polarities (Fig. nil. When the size of a bit-cell becomes 
comparable to the average diameter of grains, the effect of 
an irregular grain boundary around an edge tends to be the 
dominant source of noise that degrades the overall detection 
performance of a system. The proposed model includes the 
effect of the irregular grain boundary by introducing random 
variables corresponding to each edge in the writable area. If 
Xi and Xj have the same value (i.e., same polarity), the term 
expressing the jitter-like noise \xi — Xj\zij disappears. On 
the other hand, when Xi and Xj have opposite polarities, the 
term \xi —Xj\zij influences the received signal yi as a signal- 
dependent noise. 

For the following analysis, it is convenient to introduce a 
vector notation y = Ax -f q(x) to express the model (l8]l of 
the received signals, where 

K= {xi,X 2 ,...,Xn)^, y = {yi,y 2 ,---,ynf ■ (9) 

The matrix A G represents the linear interference 


q(x) = (( 7 i(x),q 2 (x),...,(?„(x)), (10) 

where the explicit form of qi(x) is given by 

(7i(x) = ^ \xi - Xj\zi^j + Zi. (11) 


The first and second moments of the noise term qi{x.){i G 
[l,n]) can be immediately derived as follows. The second 
moment of qi{x.) is evaluated as 

^ (xi-f (12) 

kGN(i) 

= ^ (Xi -Xfe)^cr| -f cr^. (13) 

kGN(i) 


In the above derivation, the independences of Zij and Zi 
are exploited. In a similar manner, the second cross moment 
E[q,{x)qj(pfi]{i fj) is given by 


E[q,{x)qj{x)] 


{xi — Xj)^aj, Bi and Bj are adjacent 
0 , otherwise. 


In summary, the covariance matrix regarding q(x), which is 
denoted by 5'(x) = {S'ij(x)}, is given by 




{xi-Xj)‘^a], 

0 , otherwise. 

(14) 


The vector q(x) can be obtained by multiplying a real 
matrix to a vector of independent Gaussian random variables. 
This implies that q(x) obeys the multi-dimensional Gaussian 
PDF. We therefore have the conditional PDF given by 


-PY|x(y|x) 


1 

X exp - ^x)'^S'"^(x)(y - Ax)^ 


(15) 


It is trivial to see that the conditional PDF of the two-bit-cell 
model is a special case of this conditional PDF (fTSl l. 


As an example, let us consider the case in which the writing 
area consists of 2 x 2 bit-cells. When all four bit-cells have the 
























same polarity, the corresponding covariance matrix becomes a 
diagonal matrix with the diagonal element because no jitter¬ 
like noise occurs. The worst pattern in terms of noise variance 
is the checker-board-like pattern (i.e., any pair of adjacent bit- 
cells have distinct polarity) because every edge between two 
bits causes an irregular grain boundary effect. In such a case, 
the covariance matrix has the following form: 
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where the conditional PDF fV|x(y|x) is given in (fTsT l. The 
entropy of Y is thus given by the following multi-dimensional 
integration: 

HiY) = -f PY(y)log2-PY(y)dy- (19) 
J 

The entropy iF of a multi-dimensional Gaussian random 
variable with a covariance matrix K 0 is given by 

7J = (l/2)log2((27re)'^|iT|). (20) 

Using this formula, we can derive the conditional entropy 
iF(Y|X) as follows: 


IV. Evaluation of symmetric mutual information 

Evaluation of the capacity of the simplified TDMR channel 
is of practical importance. The capacity indicates the areal 
density limit for the TDMR channel without making any 
assumptions on coding and detection methods and may also 
give us insight into the design of a read-back system including 
detection algorithms and coding. 

In the previous section, we have derived the closed form 
expression of the conditional PDE that perfectly character¬ 
izes the channel model. Therefore, in principle, we have 
sufficient information to evaluate the capacity of the channel 
given as C = maxp^ I{'X.;Y). The random variables X = 
(Xi,X 2 , ..., Xn)'^, Y = (Yi, Y 2 , ■ ■ ■, Yn)'^ represent written 
signals in a writable area and read-back signals, respectively. 
The symbol /(X; Y) is the mutual information defined by 

/(X;Y) = iJ(Y)-iJ(Y|X). (17) 

However, the problem of evaluating the capacity of this 
channel is not trivial. Of course, maximization of mutual 
information /(X; Y) in terms of the prior distribution Px is 
a computationally difficult problem. Even the evaluation of 
the mutual information itself is not a simple problem because 
we need to handle signal-dependent noise that results in a 
nonuniform conditional PDE. Another computational difficulty 
comes from the intrinsically high dimensionality of the model. 
High-dimensional integrations are required for evaluating the 
mutual information. 

In the present paper, we will focus on the symmetric mutual 
information instead of the channel capacity. The symmetric 
mutual information is the mutual information under the as¬ 
sumption that the written signal X is equiprobable. Erom the 
definition of the capacity, it is clear that the symmetric mutual 
information is a lower bound of the capacity and can be used 
as an approximate value of the capacity. This simplification 
makes the problem computationally tractable. In this section, 
we present a Monte Carlo method for evaluating the symmetric 
mutual information of the simplified TDMR channel model. 


A. Symmetric mutual information 

Under the assumption that X is equiprobable, the proba¬ 
bility Pviy) is given by 


Pviy) = Yl ^Y|x(y|x)Px(x) 
xe{+i,-i}" 

= 2-" ^ PY|x(y|x), 

xe{+i.-i}" 


H{Y\X) = -^Px(x) f PY|x(y|x)log2PY|x(y|x)dy 

X “'y 

= 2-"^ilog2((27re)-|P(x)|). (21) 

X 

The range of summation x S {-fl,—l}"is hereafter omitted. 
The symmetric mutual information, denoted by Is, is thus 
given by 


Is 


/yGR" 


Ex-fV|x(y|x) 


l0g2 


Ex^Y|x(y|x) 


- 2-"^-log2((27re)"|^(x)|). 


dy 

( 22 ) 


B. Monte Carlo method for evaluation of Is 

In order to circumvent the numerical difficulty of high¬ 
dimensional numerical integration, we use a Monte Carlo 
method to evaluate the symmetric mutual information. The 
pseudo code of the Monte Carlo method is summarized as 
Algorithm 1. 


Algorithm 1 Monte Carlo method for evaluation of Is 

I: According to (l2lT i. evaluate i7(Y|X) exactly. 

2: t := 0;s := 0 
3: while t < tjnax dO 

4: Generate x according to the uniform distribution. 

5: Generate a received vector y = Ax + q{x). 

6: According to (fTSl l. evaluate Pviy) exactly. 

7 : S := S-l 0 g 2 PY(y) 

8: t:=t + l 

9: end while 
10: 0 :— sjtjjiax 

11: Output 0 — i7(Y|X) as an estimate of Is- 


The parameter t^ax is a given positive number that indi¬ 
cates the number of iterations. In the following, we explain 
the important steps in Algorithm 1. In line 1 in Algorithm 1, 
i7(Y|X) is evaluated according to (1211 1. Since the number of 
summands of (|2TI) is 2", this step requires time complexity 
(9(2"). The step at line 6 is also a time consuming part of the 
algorithm that requires (9(2")-time for evaluation. Note that 
the evaluation of Pviy) is repeated tmax times. Thus, this 
part dominates the overall computation time. The quantity 6 
is an approximate value of T7(Y) as 

0^Ep,[-\og,P^{y)]=H{Y). 


(18) 










There is a tradeoff relationship between the computation com¬ 
plexity and accuracy of the approximation regarding 9. The 
accuracy improves as the number of iterations tmax increases 
but a larger tmax results in a longer computation time. 

V. Numerical results 

In this section, several numerical results obtained by the 
Monte Carlo method proposed in the previous section will 
be presented. Figure Q shows the plots of symmetric mutual 
information as a function of the standard deviation of the 
system noise a^. The curves of two cases such that aj = 0.4 
and tTj = 0.8 are shown in Fig. [T] In this case, numerical 
integration can be performed in order to obtain the symmetric 
mutual information because the number of bit-cells is small. 
We thus plot the results of numerical integration and the results 
of the Monte Carlo method in Fig. |7] The number of iterations 
is set to tmax = 10000 in the Monte Carlo simulations. 
Although certain statistical fluctuations are observed for the 
curves obtained by the Monte Carlo simulations, the curves 
obtained by the two methods are in reasonable agreement. We 
can also observe a tendency whereby the symmetric mutual 
information decreases as aj increases. 



Fig. 7. Symmetric mutual information of two bit-cell model as a function 
of (Js (numerical integration and Monte Carlo method) 

Figure |8] includes the curves of the symmetric mutual 
information of a simplified TDMR channel consisting of nine 
bit-cells. The interference coefficients a = 1.0,/3 = 0.5 are 
assumed. The number of iterations is set to tmax = 1000. 
In the case of aj = 0.8, the symmetric mutual information 
rate (mutual information divided by the number of bit-cells) 
approaches 2/3 at as = 0.3. 

VI. Conclusion 

In order to capture the qualitative features of the writing 
and read-back processes of TDMR systems, we proposed 
a simplified TDMR channel model. The simplicity of the 
proposed model enable us to derive the closed from of the 
conditional PDF representing the probabilistic nature of the 
channel. The conditional PDF is Gaussian distributed and 
is parameterized by the signal-dependent covariance matrix. 
Based on this conditional PDF, we developed a Monte Carlo 



Fig. 8. Symmetric mutual information of simplified TDMR channel model 
as a function of as (9 bit-cells, a = 1.0,/3 = 0.5, tmax = 1000) 

method for approximating the symmetric mutual information. 
The symmetric mutual information is closely related to the 
areal density limit for a TDMR system without making any 
assumptions on coding or detection schemes. The numerical 
results presented herein suggest that a low-rate coding, such 
as 2/3 or 1/2, may be necessary when jitter-like noise becomes 
dominant. The channel model is useful not only for estimation 
of the capacity of TDMR systems but also for the development 
of efficient detection algorithms and for the design of two- 
dimensional codes suitable for a TDMR channel. 
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